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ABSTRACT 

> 

(N 

On . Two selected regions in the molecular gas spiral arms in M51 were mapped 

Q ■ with the Owens Valley Radio Observatory (OVRO) mm- interferometer in the 

|> ' 1^00(2-1), 13CO(1-0), 01*0(1-0), HCN(l-O) andHCO+(l-0) e mission lines. 



The C O data have been combined with the ^^00(1 — 0) data from lAalto et al.l 



(1l999l ) covering the central 3.5 kpc to study the physical properties of the molec- 
ular gas. All CO data cubes were short spacing corrected using IRAM 30m 
^ ■ (i^CO(l — 0): NRO 45m) single dish data. A large velocity gradient (LVG) anal- 

j^ ■ ysis finds that the giant molecular clouds (GMCs) are similar to Galactic GMCs 

when studied at 180 pc (120 pc) resolution with an average kinetic temperature 
of Tkin = 20(16) K and H2 density of n(H2) = 120(240) cm~'^ when assuming 
virialized clouds (a constant velocity gradient ^). The associated conversion 
factor between H2 mass and CO luminosity is close to the Galactic value for 
most regions analyzed. Our findings suggest that the CMC population in the 
spiral arms of M 51 is similar to those of the Milky Way and therefore the strong 
star formation occurring in the spiral arms has no strong impact on the molecular 
gas in the spiral arms. Extinction inferred from the derived H2 column density is 
very high [Ay about 15 - 30 mag), about a factor of 5-10 higher than the average 
value derived toward HII regions. Thus a significant fraction of the ongoing star 
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formation could be hidden inside the dust lanes of the spiral arms. A comparison 
of MIPS 24/im and Ha data, however, suggests that this is not the case and 
most of the GMCs studied here are not (yet) forming stars. We also present low 
(4.5") resolution OVRO maps of the HCN(l-O) and HCO+(1-0) emission at the 
location of the brightest ^^C0(1 — 0) peak. 



Subject headings: galaxies: ISM 
galaxies: individual(NGC 5194) 



ISM: molecules — radio lines: ISM 



1. Introduction 



Since the first discovery of giant molecular clouds (GMCs) in our Galaxy, a major aim 
has been to unders tand their range of physical properties and their dependence on envi- 
ronment. Recently, iBolatto et al.l ( l2008l ) showed that the GMC properties are remarkably 
similar over a large range of environments using a large set of interferometric observations 
of local dwarf galaxies as well as M31 and M33. Most studies of spiral galaxies outside the 
Local Group were undertaken using single dish observations where the GMCs are no longer 
resolved. Thus it is not obvious if the results of these studies can be directly applied to 
the GMC population present in these galaxies. Detailed analysis of the molecular gas at 
roughly GMC-scale resolution we re only obtained for the centers of a iew nearby galaxies 
such as IC 342, Maffei 2 and M 82 jMeier fc TurneilboOll : iMeier et al.lboosi : IWeiB et al.lboOll ) 
all suggesting that the properties of these GMCs differ from those found in galaxy disks by 
exhibiting higher kinetic temperatures and lower conversion factors likely due to the effect of 
the more vigorous star formation present in the centers of these galaxies. Thus it would be 
interesting to know whether GMCs residing in galactic disks with enhanced star formation 
are similar to local GMCs or closer to those found in the central regions of external galaxies. 

The Whi rlpool galaxy M 51 is an ideal target being one of the closest (D=8.4Mpc with 
1" ~ 40.7pc; iFeldmeier et al.lll997 ) almost face-on grand-design spiral gala xies. A significant 



fract i on of its molecular gas is found in the spectacular spir al arms (e.g. iScoville fc Young 
19831 : lAalto et al.l Il999l : IShetty et al.l 120071 : iKoda et al.l 120091 ) . Therefore it offers an unique 
opportunity to study the physica l properties of the molecular gas within the spiral arms. 
For example, iKramer et al.l (120051 ) combined single dish data with observations from ISO to 
investigate the state of the interstellar medium (ISM) in the center and selected arm regions 
finding indications that the ISM in M51 might differ from that of the Galaxy as most of 
the [CII] line emission in M 51 can be explained by arising from a photon-dominated region 
(PDR) unlike the 50% assumed for the Galaxy. Due to M51's low inclination [i ~ 20°) 
non-circular motions can be more easily separated from the galactic rotation. The unusually 
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large streaming motions (60 -150 km s ^) found in M51 imply a very strong density wave 
whose general features ar e velocity discontinuities and streaming across the spiral arms (e.g. 
Roberts fc Steward EqStI i 



Besides the large number of Giant M olecular Cloud Associations (GMAs) (up to 16 
within one spiral arm, lAalto et al.l 119991 ). numerous star clusters an d HIT regions resid e 
in the spiral arms as is obvious in high resolution HST imaging (e.g. IScoville et al.l 120011) . 



Bastian et al. 



2005 



These populations have been extensively studied using HST data (e.g. 
Lee et al.ll2005l : IScheepmaker et al.ll2009l and references therein). Two interesting results in 
the context of this paper are that younger clusters show higher extinction and are located 
closer to spiral arms. 

The paper is organized as follows, the observations are presented in §2] while the data 
are described in §21 A non-LTE analysis of the spectral lines is given in §H The relation 
between the CMC properties and star formation is investigated in ^ A general discussion 
(^ and a summary of the findings {^ are provided at the end. 



2. Observations 



2.1. Millimeter Interferometric Observations 



For this multi-transition study of the molecular gas spiral arms in M 51a, two prominent 
regions (Fig. [T]) were selected based on the interferometric map of the ^^C0(1 — 0) line emis- 
sion obtained earli er with the Owens Valley Radio Observatory (OVRO) mm-interferometer 
( lAalto et al.lll999l ): (a) the region containing the brightest ^^C0(1 — 0) peak (hereafter: 
west), and (b) the region showing evidence for shear across a spiral arm south of the nu- 
cleus (hereafter: south) comprised of the GMAs A8, A9, and AlO (nomenclature following 
Rand fc Kulkarnilll990l ). Each region was covered by a single pointing with the OVRO inter- 
ferometer, and both regions are in the Ml spiral arm. The coordinates of the pointing centers 
are given in Tab. [1] We observed the following key tracers for molecular gas: ^^C0(2 — 1), 
^^CO(l-O), 0^^0(1-0), HCN(l-O) and HCO+(1-0). In order to accommodate simultaneous 
observations of the lines of 1^00(1-0) and 0^^0(1-0) as well as HCN(l-O) and HCO+(1-0), 
the spectral set-up for the southern and western pointing was centered at slightly different 
velocities. Tab. [2] lists the velocity of the central channel as well as the resolution adopted 
for the data analysis. The data were obtained between January 2001 and March 2002 in 
the C, L, E, and H configuration (for details see Tab. E]). J1159+292 and/or J1310-t-323 
served as phase calibrators and were observed every 20min between source observations. For 
passband calibration 3C273, 3C84, and/or 3C345 were used. Observations of Neptune or 
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Uranus were observed in several tracks throughout the observing period and allowed the flux 
calibration of the phase calibrator which was when interpolated for the M51 observations. 
The uncertainty of the absolute flux calibration is about 10% (15%) at Srnni (Imni). Th e 



calibration of the uv data was done using the OVRO software MM A (IScoville et al.l 1994I ). 



For mapping and CLEANing we used the software package MIRIAD (jSault et al.lll995l ). The 



FWHM of the primary beam of the OVRO interferometer is about 60" (38") at 3mm (1mm). 
In add i tion t o the new observations, we use the ^^C0(1 — 0) mosaic obtained by 



Aalto et al.l ( 1l999l ). The mosaicked map consists of 19 pointings observed in the C, L and 
H configuration with OVRO. The final spectral resolution is 7.8 kms~^ with 33 channels 
centered at a velocity of vlsr = 472 kms~^. The sensitivity of the naturally (uniform) 
weighted data is 19mJybeam"^ (25mJybeam~^ ) with a beam of 3 .9" x 3.3" (2.9" x 2.1"). 



A complete description of the dataset is given by lAalto et al.l (119991 ) 



2.2. Millimeter Single Dish Observations 

For a meaningful comparison of the different tracers of the molecular gas, the interfero- 
metric data need to be corrected for missing short spacings (SSC - short spacing correction). 
Therefore, we observed the regions of the western and southern OVRO pointings with the 
IRAM 30m to obtain total power data. Around each pointing center a region of 80" x 80" 
was mapped in the On-The-Fly (OTF) mode in July 2001. The observations were carried 
out using the AB receiver combination in two frequency setups covering simultaneously the 
1^00(1 - 0) and ^300(2 - 1) or the 0^^0(1 - 0) and i3C0(2 - 1) lines in dual-polarization 
mode at 3 and 1mm. Since we did not observe the ^^C0(2 — 1) line at OVRO, the 1mm 
data will not be described in the following. 

The telescope beam size at the 3mm observing frequencies is 22". We used four of the 
1 MHz filter-banks as spectral back-ends. Each unit has 256 channels resulting in a velocity 
resolution of 2.6 kms~^ and a total velocity coverage of 670 kms~^ at 3mm, respectively. 
System temperatures were typically ^170K (T^) for both frequency setups. Pointing was 
checked frequently and was found to be stable to within 3". Calibration was done every 
15 min using the hot/cold-load absorber measurements. The OTF observations were done 
in equatorial coordinates, alternating the two orthogonal scanning directions. The scan 
speed was 2"/s with a dump time of 2 s and steps orthogonal to the scanning direction of 4" 
yielding spectra spaced by 4" on sky. Data were reduced using the CLASS software. Linear 
baselines were subtracted from each spectrum and intensities were converted to main beam 
brightness temperatures (Tmb)- The calibrated, baseline-subtracted spectra were used to 
produce a 3-D data cube. During the gridding we reduced the velocity resolution to 5kms~^ 
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for both lines. The resulting rms noise level for ^^C0(1 — 0) is 9mK in both regions. For 
the C^®0(1 — 0) line we achieved an rms noise of 6niK in the west and 12 niK in the south. 
We estimate fluxes are accurate to ±10% . 



For the SSC of the OVRO mosaic of the i2C0(l - 0) data from lAalto et all fll999h and 
our OVRO 1^00(2 - 1) observations, we used the ^^00(1 - 0) a nd 1^00(2 - 1) single dish 
maps obtained with the NRO 45 ni telescope (jNakai et al.l Il994[ l^l and the HERA receiver 
array at the IRAM 30m telescope (ISchuster et al.l 120071 ). respectively. The ^^C0(1 — 0) data 
cube has spatial resolutions of 16" and a spatial gridding of 7.5" on sky for the region of 
interest here. The ^^C0(2 — 1) data cube was observed in OTF (On-The-Fly) mode and 
has a spatial resolution of 11" and a spatial gridding of 7". The rms noise levels at the 
resampled spectral resolution of 5.3kms~^ and 4.5kms~^ is 60mK and 30 mK (Tmb) for the 
^^C0(1 — 0) and ^^C0(2 — 1) line, respectively. We have checked and corrected the astrometry 
of the single-dish maps relative to the OVRO maps by comparing the peaks in the integrated 
intensity maps before the short spacing correction. For all maps the astrometry turned out 
to be consistent within the pointing errors of the single-dish observations with applied shifts 
of {5RA = 4.5",5DEC = -1.5") and {6RA = 3.2",5DEC = 0.0") for the i2C0(l - 0) and 
^^C0(2 — 1) line, respectively. 



2.3. Short Spacing Correction 

In order to recover the flux resolved out by the interferometer (missing flux), we com- 
bined the OVRO interfero metric and s ingle dish data following the short spacing correction 
(SSC) method outlined by lWeiB et al.l (120011 ). In brief, this method uses the cleaned interfer- 
ometer map as a starting point and replaces the central part of the uv-plane with visibilities 
calculated from the single dish map using the same spatial and velocity grid and flux units 
as the interferometer map. The method requires knowledge of the interferometer clean beam 
and the single dish beam (which has been approximated by a Gaussian). The method has 
no free parameters except the choice which part of the u^^plane is replaced by the single 
dish data. All visibilities shorter than the smallest projected baseline of the interferometer 
data were taken from Fourier transforms of the single dish data cubes. For the analysis, we 
produced short spacing corrected data cubes at 4.5" (2.9") for both (west only) regions. We 
flnd that the missing flux in the OVRO pointings is ~ 50, 40 and 20% for the ^^C0(1 — 0), 
1^00(1 - 0) and 0^80(1 - 0) lines, respectively, and ~ 70% for the 1^00(2 - 1) line. The 
missing flux percentage is lower towards the strong emission peaks in the interferometer 



^ http : //www . nro . nao . ac . jp/~nro45inrt/C0atlas/ 
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maps. All maps presented in the following are the short spacing corrected versions. 



2.4. HST Archival Data 

For comparison to sites of recent star formation we use archival HST data from the 
W FPC2 camera i n the V band and Ha line as well as the NICMOS Paa line as presented 



bv lScoville et all (120011 ) 



3. The molecular ISM 

In order to study the molecular gas emission, we created moment maps from the short 
spacing corrected data cubes using the GIPSY (Groningen Image Processing System) task 
'moment'. We used a Sex (5cr) or 5cr {7a) limit and required that emission above this threshold 
is present in at least two adjacent channels to minimize the noise for the intensity maps 
(velocity fields). Note that this is a very strict cut and low level emission will be missed by 
the moment maps presented here. 



3.1. CO Morphology 

Strong line emission has been detected from the ^^C0(1— 0), ^^C0(2— 1), and ^^C0(1— 0) 
transition while 0^*^0(1 —0) is only detected at a low level in some smaller peaks (see Fig. [2l 
[3]and|l]). However, the relative brightness of peaks within the arms is different for each line 
suggesting that different physical conditions might be present in different cloud structures. 
A slight change in CO morphology is apparent for the western region when comparing the 
i2C0(l - 0), i2C0(2 - 1), and i3C0(l - 0) line emission at 4.5" (Fig. ED and 2.9" (Fig. 
[3]) resolution: the small-scale geometry like the ridge connecting the peaks at a relative 
declination at +5" and 0" becomes more evident in the ^^C0(1 — 0) and ^'^C0(1 — 0) maps. 
The difference is likely due to an extended low excitation component which is contributing 
more to the ^^C0(1 — 0) line emission than the ^^C0(2 — 1) one. 

No significant offset is seen between the peaks in ^^C0(1 — 0) and ^'^C0(1 — 0) (Fig. [5] 
andEj). We verified that the location of the peaks in the moment maps is not significantly 
affected by the inclusion of the single dish data. Small offsets of a fraction of the beam 
are seen for some peaks (mainly in the ^^C0(1 — 0) moment maps) as is expected in the 
case of low S/N data (being boosted by the additional signal from the single dish data). In 
addition, as the physical conditions across the spiral arms are changing as shown in §11 a 
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perfect one-to-one correspondence between the ^^C0(1 — 0) and ^^C0(1 — 0) morphology 
is not expected . The absence of an obvious shift is in apparent contradiction to the result 
of IXosaki et aD J2002I ) who found that the ^^00(1 - 0) emission is leading the ^^00(1 - 0) 
emission in particular in the southern region. One explanation could be that this apparent 
offset is caused more by an offset in the diffuse gas than in the more clumpy components 
which are better resolved with the OVRO data. Another possibility is the instrumental 
offset we had to correct for when aligning the single dish and interferometric data (see §2.2p . 
As the offset is mainly in RA for the ^^C0(1 — 0) single dish data this could explain the 
observed difference seen for the Western arm. If an instrumental effect is also causing the 
mismatch for the Southern arm, the astrometry of the ^^C0(1 — 0) single dish data should 
be responsible. 



The smooth ridges present in the ^^C0(1 — 0) intensity map of lAalto et al.l (119991 ) break 
up into smaller peaks at ~1.6" resolution achieved in the ^^C0(2 — 1) data. In particular the 
western region (Fig. [7]) splits into two gas lanes over a length of about 10" (~ 410 pc) while 
the southern region (Fig. |8]) show s up to three s purs a long the leading side of the spiral arm 
similar to the features studied by ICorder et al.l (120081 ). The peaks in both regions appear to 
be slightly resolved at that resolution. 



3.2. CO Kinematics 



The 1.6" resolution ^^C0(2 — 1) data allows us to study the molecular gas kinematics 
at about 65 pc resolution. Alth ough the overall k inematics corresponds well with the low 
angular resolution velocities of (jAalto et al.l Il999l . their Fig. 3), the iso- velocities showing 
small deviations indicate that some GMC associations (GMAs) move faster or slower (fig. 
13 middle). While the general deviations from a pure rotational disk are present as well 
indicating an overall non-circular component, some peaks stand out in the velocity dispersion 
map suggesting that a few peculiar velocities in the velocity maps are due to a superposition 
of GMCs rather than actual streaming motions. The most obvious example is the peculiar 



veloci ty field structure identified as streaming motion across the southern arm by lAalto et al. 



(119991 ) that is likely caused by a superposition of two GMCs rather than a distinct velocity 
feature due streaming (Fig. |8] middle and bottom). The high velocity dispersion observed in 



the G MA at offsets of +5" and +0" (GMA AlO following the nomenclature of lAalto et al. 



19991 ) and the break-up into several sub- components of the GMA at offsets +12" and +2" 



(GMA A9) support our interpretation. In addition, the 1.6" resolution ^^C0( 2 — 1) data 
reveal s that the higher velocity dispersion ridge seen in the ^^C0(1 — 0) mosaic of lAalto et al. 



( 119991 ) in particular in the western arm belongs to the second set of line peaks on the inner 



side of the spiral arm (see Fig. [7] right). 



3.3. Dense Gas Tracers 

HCN(l-O) and HCO+(1-0) hne emission was only detected in the western arm region at 
a low significance level (see Fig. |2]). HCO+(1-0) is confined to the bright CO peak at about 
-2" and +0" while HCN(l-O) might be associated also with another bright CO peak north 
of this position. Given the low significance of the detection these data are not included in 
the following analysis. 



4. Non-LTE Analysis 

In order to estimate the physical conditions in the molecular gas along the gas spiral arm 
we measured the CO line emission in all four transitions for several positions along the spiral 
arms. As emission from the C^^0(1 — 0) line was only detected in a small fraction of the 
positions, we will focus on the ^^C0(1 —0), ^^C0(2 — 1), and ^^C0(1 — 0) line measurements 
for our analysis. 



4.1. Line Ratios 

14 and 9 positions were selected in the western arm region by basically centering them 
on distinct peaks of the CO line emission in the 2.9" and 4.5" resolution maps, respectively. 
In the southern region 7 positions were selected in a similar manner at 4.5" resolution only. 
The western 2.9" resolution and southern 4.5" resolution positions are indicated in Fig. M 
and [TOl respectively. For the measurement of the line ratios R21 of the ^^C0(2 — 1) to 
i2C0(l - 0) line, Ri3 of the i2C0(l - 0) and 1^00(1 - 0) line and Ris of the i2C0(l - 0) 
and C^^0(1 — 0) line all data cubes were converted from flux density with units of Jy/beam 
to brightness temperature with units of K. For each position spectra of the four CO lines 
were extracted at the central pixel. In a few positions spectra showed evidence for two line 
components suggesting that more than one CMC is contributing to the emission peak. In 
these cases two Gaussians were fit to the spectra. 

The integrated line fluxes were derived from Gaussian fits to the spectra and are listed 
in Tab. [5]and Tab. [6]for the 2.9" and 4.5" resolution data. Upper limits for the line flux were 
derived using three times the rms in the spectrum multiplied by the FWHM measured for 
the ^2C0(1 -0) line. The listed line width corresponds to the FWHM of the ^^00(1 -0) line. 



-9- 



The line ratios i?2i, -R13 and -Ris were derived by scaling the overall shape of the ^^C0(2 — 1), 
^•^00(1 — 0) and 0^*^0(1 — 0) spectra to the corresponding ^^C0(1 — 0) spectrum and are 
listed in Tab. O and M as well. Thus the listed line ratios might not agree perfectly with a 
line ratio derived from the integrated line fluxes. This approach is more robust in the case of 
low S/N spectra than using the peak or integrated line flux which could be strongly effected 
by noise peaks. If no hne was detectable in the data, we used the peak flux (upper limit) of 
the lines in question to derive a lower limit for the R13 and Ris ratios. 

The derived line ratios R21 and R13 for the 2.9" resolution western and 4.5" resolution 
southern region are shown as a function of position in Fig. [9] and [TOl respectively. In the 
western region the R21 ratio is basically monotonically decreasing as a function of galacto- 
centric distance with a value of i?2i=0.72 at position wl corresponding to a radius of about 
1.1 kpc (~ 27") and values between 0.42 and 0.47 for positions wll to wl4 which are at a 
radial distance of about 1.5 kpc (~ 37"). The -R13 ratio on the other hand shows a spread 
of values between 4.5 (wl4) and 7.5 (w5) with no clear radial trend. However, most of the 
higher values are found close to the area of brightest CO emission. While the western region 
covers a larger range in galactic radii, the southern region is basically tracing molecular 
gas located at roughly 2.6 kpc distance from the center. Thus the R21 and R13 have fairly 
similar values over most of the positions, except for the most eastern ones (S6 and S7) which 
are located close to a very large HII region (see ^ and show a double line profile. It is 
interesting to note that both regions exhibit one component where the ratios are slightly 
higher than those found for the southern most positions in the western region. We note that 
all trends seen within a single pointing should be real, as calibration uncertainties should 
only affect the absolute flux levels and thus all positions investigated in the same way. 



4.2. LVG Parameters 

To constrain the physical conditions of the GMCs we use an spherical, isothermal one- 



comp onent large velocity gradient (LVG) model (IScoville &: Solomonlll974l : lGoldreich fc Kwan 



19741 ). This treatment of the radiative transfer is valid for clouds in which a large-scale 
velocity gradient effectively makes the radiative transfer a local problem not requiring self- 
consistent solutions for the entire cloud. The assumption for this model is that due to 
the velocity gradient present the emission of a molecule is Doppler-shifted with respect to 
molecules in the rest of the cloud. For the models presented here we fixed the ^^CO abun- 
dance relative to H2 and the ^^CO/^^CO abundance ratio to 8.0e-5 and 30, respectively. The 
numbers correspond to values found in GMCs in th e inner 3 kpc of the disk of the Milky 



Way (e.g. iLanger &: Penziaslll993l : iMilam et al.ll2005l ). LVG models were calculated for two 
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different assumption on the velocity gradient: 

typically for local clouds and reasonable for M51's cloudcl, b) assuming a virialized cloud 



a fixed value of ^ set to Ikms ^pc ^ 



such that ^ 

dr 



3.1 



n{H2) 
104 



[e.g. 



Goldsmith! 120011 ) . LVG line intensity ratios were computed 



for kinetic gas temperatures Tkjn betwee n 3 and 200 K an d H2 densities logn(H2) between 
2.0 and 4.0 (for details of the model, see lWeifiet aPboOlh . 



4.3. M51 model results 

The results of the LVG analysis are presented in Fig. IH] and [10] and summarized in Tab. 
[71 Typical relative errors on the derived temperature are of the order of (75-100)% implying 
that solutions preferring higher temperatures are in principle also consistent with the low 
temperatures found for several locations. For the H2 densities of A^(H2) and n{B.2) we find 
typical relative errors of about 25% (going up to ~50% in the less constrained cases). These 
errors are independent on the assumptions made for the velocity gradient. 

The bulk of the CO emission probed here arises from cold clouds (Tkin ~ 16-20 K, 
depending on the assumption made for the velocity gradient) with moderate H2 density of 
n(H2) ~ 240 — 120 cm~^. It is interesting to note that the H2 density in the assumption 
of a virialized cloud is typically (50-65)% of the H2 density derived for a fixed velocity 
gradient. The average conversion factor derived from the LVG analysis of Xco = (1.3—2.0) x 



CO. 



:i.8±0.3 ) xlO^o cm-^K-^km-^E 



10^° cm~^K^^km^^s is close to the Galactic valueoof X, 
This value was recently derived for the Milky Way by iDame et al.l (l200l[ ) suggesting that 
the GMCs analyzed here are similar to those in the Milky Way, especially when assuming 
that they are all virialized clouds. Thus the data presented here mainly probe the diffuse 
component of the ISM and not the dense gas that is involved in current star formation. In 
the remainder we refer to the LVG results assuming virialized clouds. However, all trends 
discussed are present in both scenarios. 

The most obvious spatial variation of the gas physical condition is a decrease of the 
kinetic temperature with increasing galacto-centric radius along the western spiral arm re- 
gion. Temperatures decrease from about 27 K at 1.1 kpc (position wl) to <10K at 1.6 kpc 
(position wl4). This might indicate an effect of the interstellar radiation field (ISRF) on the 
cloud properties as these radii probe the transition from the central (bulge) region into the 



^For an average line width of 40knis ^ (see e.g. Tab. [5]) and an average cloud complex size of ~50pc 
(see ©. 

•^First determined bv lSanders et al.l (|l984l) . 
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disk (for a discussion of the effect of HIT regions see ||5]). In the southern region, with an 
almost constant galacto-centric radius of ~ 2.6 kpc, no such gradient is detected. These more 
distant regions are on average 5-10 K warmer than the furthest regi ons in the Weste r n arm 



(position wl2 and wl4). The derived radial temperature profile by iMeijerink et al.l (|2005[ ) 
used to model the heating mechanism at 850/im shows a continuous decline of temperature 
with radius, howeve r, the change over the dist ances relevant here is less than 2 K. The dust 



temperature map of iBenford fc StaguhnI (120081 ) derived by fitting a grey-body thermal emis 



sion spectrum to data at 70, 160, 350 and 850 /im shows also an enhanced dust temperature 
at the center with a decline for larger radii. However, the dust temperature in the spiral 
arms themselves is higher than in th e inter arm regions. Overall the change in temperature 



is not dramatic. iBenford fc StaguhnI (120081 ) report a mean temperature of (25±3) K for the 
disk while the inner 3 kpc region has a higher temperature of 31 K. Therefore it appears that 
more local changes (i.e. within the spiral arm) of the ISRF than the global radial trend could 
be causing the apparent observed trend in the Western arm. However, higher resolution data 
such as observations from the Herschel Space Observatory and subsequent detailed modeling 
are required to address this in more detail. 

There is no obvious large difference in the temperature and H2 density of GMCs close 
to or around star forming regions as traced by the presence of Ha emission (see §5], also 
indicated in Tab. [7]). This finding suggests that the CO emission arising from or close to 
star-forming regions is still dominated by cold gas that is surrounding the HII regions and 
not immediately affected by the enhanced radiation from the HII region. There might be 
a slight effect present as half of the positions in the Southern arm encompass HII regions 
(traced by Ha emission) and show slightly higher temperatures. The impact of star formation 
on the molecular clouds (e.g. heating at the cloud surfaces) is not measurable at a linear 
resolution of 180 pc (or even for the higher 120 pc resolution data). Only two positions stand 
clearly out: w5 in the western region and S7(B) in the southern region that show the highest 
kinetic temperature of T^in = 50 K, average H2 densities and the lowest conversion factor 
of Xco = 0.5(0.8) X 10^° cm~^K~^km~^s. Both GMC complexes are located on the inner 
edge of the CO spiral arm suggesting that they either live in a different environment than 
the clouds studied in the other locations or that their ISM has been heavily affected by the 
ongoing neighboring massive star formation [in the case of S7(B)]. 



4.4. Xco measurements for M51 

Several studies were conducted in the past to derive the conversion factor Xco for the 
molecular gas in M51. Most of these studies did find values well below a Galactic value of 
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Xco = (1-8 =t 0.3) X 10^° cni~^K~^kni~^s f lDame et al.ll200ll ). Several reasons could explain 
the apparent mismatch between our results and these earlier works. One main advantage of 
the work presented here is its superior spatial resolution compared to the previous results. 



Usi ng single dish observ a tions of ^^CO and ^^CO line transitions and three different 
methods ICarcia-Burillo et all Jl993h find a value of Xco = (0.8 ±0.3) x lO^Ocm-^K-^m-^s 
for the center and the molecular spiral arms from LVG modeling and an even lower value when 
using only the ^^CO emission line data. Assuming a const ant ratio between th e extinction and 
the H2 column density N{H2) as found for our Galaxy ( Bohlin et al. 19781) gave a G alactic 
value of Xco ~ 1-8 x 10^°cm"^K~^km~^s. Extending this work iGuelin et al.l ( 119951 ) added 
1mm continuum measurements to estimate Xqo via the gas-to-dust ratio, again getting 
factors about 4 times lower than the standard value. The most likely explanation for the 
mismatch is that in these studies individual GMC complexes are not resolved and thus a 
significant mixing of vastly different cloud properties can happen. 



Using interferometric observations of the molecular gas disk iNakai fc Kund ( 1l995l ) ap- 
plied the Ay methoqj to tho se HIT regions in the spir al arms that had both Ha as well 
as radio continuum emission (Ivan der Hulst et al.l 119881 ) that allowed for a good extinction 
estimate. They obtained a value of Xco ~ (0.9 ± 0.1) x 10^°cm~^K~^km~^s at about 10" 
resolution. Among possible sources of error they list a bias for probing gas with higher exci- 
tation temperatures as they are looking toward HII region where the gas could be potentially 
more heated. They also caution that the conversion factor in the inter-arms might be com- 
pletely different as the conditions for the molecular gas are expected to vary as well. Given 
that most of the gas for which we derived the conversion factor is not close to HII regions 
(see ^ those measurements are very likely affected by probing a very different environment 
than we did. 

I nterestingly, applying vi rial r nass measur e ment s to the GMAs in the spiral arms of 
M51, iRand fc Kulkarnil (Il990[ ) and lAdler et al.l ( 119921 ) come to very opposite conclusions of 

^km~^s, respectively. Both 



X, 



CO 



3 X 10^° cm-^K-^m-^s and Xco 



1.2 X lO^Ocm-^K" 



studies measure virial masses assuming that the GMAs are resolved with typical sizes of 300- 
400 pc. In order to apply this method the assumption must be that the GMAs are bound 
complexes which might not be true given the fact that it is unclear whether they are just a gas 
accurnulation due to orb it crowding or indeed self-gravitating entities (e.g. JRand fc Kulkarni 
Adler et al.l 119921 ) . This is still an open issue discussed by iKoda et al.l ( l2009l ) for 



1990 



the GMAs identified in their high fidelity map as the elongated structure could easily be 



*The adopted method bv lNakai fc Kund ( 19951) takes into account gas located behind HII regions on a 
statisticaUy basis. 
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interpreted as being due to orbit crowding rather than resembhng a bound entity even at 
their 4" resolution. As our high resolution ^^C0(2 — 1) data shows that these GMAs break 
up into several clumps that might no longer even correspond to a single gas lanes (as e.g. 
seen in the Western region), it seems sensible that the assumption of bound systems is not 
valid for all GMAs. In addition, enhanced velocity dispersion due to shocks or streaming 
motions will critically affect the measure ments of the line w idth at the 10" resolution of the 



data used by the previous studies. Since lAdler et al.l (119921 ) mainly analyzed GMAs located 



in the central region where the dynamical properties are changing faster, this could be a 
possible explanation for the different results of both works as the method used is basically 
the same. 

Taking a different approach and c ombining measureme nts of the first three transitions of 
^^CO with those of ISM cooling lines, Kramer et al.l ( 120051 ) used a PDR (photon-dominated 



region) model. However, their derived densities for the center and two regions covering 
the spiral arms at slightly larger radii are about 1 order of magnitude higher while their 
temperatures of 12.5 and 15 K are similar to our finding s . A r ecent similar study of the 



central 3.5 kpc using a refined PDR model by iBell et al.l (J2007| ) finds a Xco of 1/10 the 
Galactic value. As both studies were also done at low resolution, it is not clear that the 
ISM cooling lines as well as the CO lines have the same distribution and filling factors which 
could affect the estimation of Xco- Given that the molecular gas studied here is clearly not 
located closely to HII regions, it is it not obvious that the PDR assumption is applicable for 
our regions (see ^. 

Since it is expected that Xco is changing as a function of environment, the assumption of 
a single excitation and a simple geometry over a large region can frequently yield misleading 
results. As our measurements have the highest spatial resolution and are confined to mainly 
probing a single environment, they should present a good estimate oi Xco for the spiral arms 
probed. Given the large uncertainties involved when averaging over large areas in a galaxies, 
they are consistent with the previous results and only highlight the fact that deriving an 
average Xco for an entire galaxy is not straight-forward. In particular the relation between 
the conversion factor and the gas density and temperature Xco ~ T]^j^n(H2)°'^ can explain 
why the conversion factor is lower for studies that include a large fraction of the (presumably) 
less dense inter-arm gas. 



4.5. Comparison to GMCs in other spiral galaxies 

The LVG models suggest that the physical conditions in the GMCs in the spiral arms 
of M51 are very similar to those in the Milky Way which have typical gas temperatures 
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and densities of 10 K and above 200 cm ^, respectively (e.g. IScoville fc Sanders! Il987l ) (for a 



recent compilation see lTielensll2005l ) . Also the observed average H2 column density in M51's 

.2) X 10^2 cm- 
Galactic plane of N(H2) ~ 4xl022cm-2 fISolomon et al.lll979h . 



GMCs of N(H2) = (2.3 - 3.2) x 10^2 cm- ^ is similar to that ty pically found for GMCs in the 



Studies of GMCs in local group galaxies found that the cloud proper t ies do vary be- 
tween clouds with and without the presence of star formation. IWilson et al.l ( 119971 ) analyzed 
seven GMCs in M33 at 22" resolution (corresponding to roughly 90 pc) using multi-line CO 
observations in conjunction with a LVG model. They find that GMCs without star for- 
mation are typically colder with Tkin ~ 10 — 20 K than GMCs with star formation present 
showing Tkin ~ 30 — 100 K for our assumed [^^cOj/f^^CO] abundance ratio of 30. While the 
temperatures for the cold CMC population is similar to our values for M51, their derived 
densities are about 1 order of magnitude higher than our estimates for M51. Interestingly, 
they note that the sphere of influence for HII regions seems limited as GMCs within 120 pc 
distance of an HII region do already show normal prop erties and no effec t of he ating. A 
recent study of GMCs in the Large Magellanic Cloud by iMinamidani et al.l ( l2008l ) also find 
a correlation between the star formation activity (measured by the Ha flux) and the tem- 
perature and density of the GMCs themselves. They propose an evolutionary sequence from 
cool (Tkin ~ 10 - 30 K) and low density (log(n(H2)) < 3) GMCs, through warm (Tkin ~ 30 - 
200 K) and cool low density GMCs to warm and dense GMCs where the first two types are 
in a young star formation phase where density has not yet reached high enough values to 
cause active massive star formation. GMCs in the last category are in a later star formation 
phase where the average density is higher. 

Thus we conclude that the bulk of the GMCs observed in M51's spiral arms are very 
similar to the ones in our own Galaxy and given their low temperatures are likely not 
impacted by (ongoing) star formation. This is not that surprising as most HII regions are 
located close to CO peak s but rarely do coincide (see ^ and as suggested by the M33 
results from IWilson et al.l ( 119971 ) the sphere of influence of an HII region is small. This 
could also suggest that GMCs located in spiral arms are di fferent from GMCs located in the 
center of starburst galaxies such as M82. IWeiB et al.l (l200l[ ) found significantly higher kinetic 
temperatures of Tkin between 50 and 190 K as well as gas densities (n(H2) between 5 x 10^ and 
2xl0^cm-^) and concluded that the ongoing massive star formation is significantly altering 
the cloud properties. 
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The Relation of the Molecular ISM to Sites of Star Formation 



In order to compare the distribution of the molecular gas w ith the sites o f recen t star 
formation we use our CO data and archival HST images from IScoville et al.l (l200l[ ). As 
expected the molecular gas coincides very well with the prominent dust lanes seen in the 
individual spiral arms. Using the 1.6" resolution ^^C0(2 — 1) data one can see that the 
emission peaks do fall into regions of high extinction obvious in the HST V band image (Fig. 
[TT] and [T2|) . In particular, the double gas lane in the brightest part of the western region 
does show a line of enhanced V band emission between the two CO lanes. 

While enhanced extinction is present across the full width of the gas spiral arms of ~ 
7" it is not uniform but rather patchy. The widths of the dust lanes visible in the HST 
V band image range from about 0.5" to 2" suggesting that the GMCs we resolve in the 
^^C0(2 — 1) data are causing the extinction and that the molecular gas is mainly distributed 
within GMCs. An upper limit on the scale height of the gas disk is about 20-80 pc if we 
assume that the extincting material can not have a larger vertical extent than a dust lane 
is wide. Combining these numbers with the derived H2 column densities N(H2) in §1] (see 
Tab. [7]) returns approximately the H2 volume densities n(H2) obtained by the LVG analysis. 
This implies that most of the GMCs seen at 1.6" resolution are mostly single entities and 
not blends of several GMCs. 

Most of the HII regions as traced by their Ha and Paa line emission lie downstream 
of the molecular gas as it is expected if star formation is occurring inside the molecular gas 
spiral arms (see Fig. [TT] and [T2l) . Since the southern region is closer to corotation (where 
the pattern speed is closer to the angular velocity) than the western region, this can explain 
why two of the ^^C0(2 — 1) peaks coincide with star formation whereas only one does so in 
the western region. Except for one case, all Paa peaks have Ha emission associated with 
them suggesting that extinction is not extreme towards HII regions. 

It is interesting to compare the average extinction to ward HII regions d erived from 
the recombination line ratios of Ha and Paa of Ay ~ 3 .0 (IScoville et al.l 120011) with those 
estimated from the CO data. We use equation (6) from iNakai fc Kund ( 19951 ) of N{H2) = 
1.87 X 10^^ Ay [mag] — N(HI)[cm~'^]/2 and our average H2 column density of < N{H2) >= 
(2.3 — 3.3) X 1 ^^ cm ~^. Since the ISM is mostly molecular in the regions studied here, e.g. 



Schuster et al.l (120071 ) find values for the ratio of surface density of J2hi I ^h^ ^ ^'^ ^'-'^ 
r< 1.6 kpc and ~ 0.5 for the southern region at 2.6 kpc, we neglect the contribution from 
atomic hydrogen. The derived extinction is about Ay ~ 12 — 18, significantly larger than 
that measured in the HII regions. This value is e ven larger (Ay ~ (15 — 33)) then using a 
Galactic relation of N {H2) / Ay = 10^^ cm~'^ K~^ s (JBohlin et al.lll978l ). Thus a lot of massive 
star formation could be hidden inside the gas spiral arms and only be detectable at longer 
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wavelengths such as the infrared and radio. Us ing the HiRes deconvo lved 24 /xm MIPS image 
(Dumas et al., subm.) and the Ha map from ICalzetti et aLl ( 120051 ) we tested this scenario 
at 2" resolution. As the 24/im emission coincides with the CO emission and is thus mainly 
arising inside (or upstream) the Ha emission (see Fig. [TT] and [T^ . this supports the picture 
that hidden star formation does occur inside the spiral arms. Using the Ha emission located 
in the spiral arms and 24/im emission arising from the same area, star formation rates (SFRs) 
can be estimateqj. Co rrecting the observe d Ha flux for the average extinction of Ay ^ 3 



found for HII region s (IScoville etal.l 120011) and using the SFR prescriptions of iKennicutt 
fll998[ ) for Ha and of iRieke et aD J2OO9I ) for the 24/im data, we find that the SFR hidden in 
the spiral arms is about 50% of that traced by Ha. This 'lack' of (embedded) star formation 
inside the gas spiral arms is consistent with the fact that the molecular gas is still cold and 
might therefore have not yet started to form new stars. 

There is no strong correlation between the location of bright HII regions and massive 
CO clouds or holes, though most HII regions tend to lie next to a ^^C0(2 — 1) peak. This 
would be in line with the picture described by iBastian et al.l (120051 ) who found a correlation 
between the radius and mass for clusters of star clusters that is similar to the relation seen for 
GMCs. They interpret this as evidence that a CMC does form an assemble of star clusters. 



6. The impact of the star formation on the ISM 

The physical properties of the CMC complexes located in the spiral arms of M51 are 
very similar to those observed for GMCs in our Galaxy. Our LVG analysis suggests that the 
conversion factor between CO line emission and H2 mass is similar to the standard conversion 
factor Xqo- This result is in apparent contradiction to other studies conducted in M51 at 
significantly lower angular resolution. The simplest explanation is that the conversion factor 
is not constant across the disk (or even within the spiral arms) but a function of environment 
(even for the diffuse gas probed here). Therefore a spatial resolution matched to the size 
of CMC complexes is required to obtain reasonable numbers. Further, the assumption of a 
single value is not correct when measuring the total amount of gas present in a nearby spiral 
galaxy. Our 65 pc resolution ^^C0(2 — 1) data shows that the spiral arms break up into 
clumps at these scales that are roughly aligned along the spiral arms. In the western region 
the CO peaks are more abundant toward the leading side of the spiral showing a slight shift 
from the peak of the more smooth emission as traced by, e.g., the 4.5" resolution ^^C0(1 —0) 



^We do not subtract the contribution from the diffuse emission as star formation appears to be the 
dominant source for the Ha and 24/im emission in the spiral arms. 
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emission (see Fig. [2] and [7]). No offset is obvious for the southern region (see Fig. [Hand [8]) 
as expected as the pattern is moving close to the angular velocity at these radii. 

At our 120 pc resolution we find no evidence that star formation in HII regions is impact- 
ing the properties of the GMCs t hat a re at about >1" (40 pc) distance, in agreement with 



findings in M33 by IWilson et al.l (119971 ). Only two positions which are located at the inner 
edge of the spiral arm exhibit elevated temperatures that may be related to the large-scale 
dynamics. 

Based on the observations presented here, we conclude that most of the GMCs in the 
spiral arms of M51 are cold, dense structures that are not heavily affected by the neighboring 
star forming HII regions. Despite the inferred high extinction of Ay ~ (15 — 30) towards 
these GMCs only about 50% of the star formation traced by Ha emission is coinciding with 
the gas spiral arms indicating that most of the molecular gas seen is not yet actively forming 
stars. 



7. Summary and Conclusions 

M51 offers the ideal environment to study the physical properties of the molecular 
gas in the spiral arms of a disk galaxy. Combining short spacing corrected interferometric 
observations of the ^^C0(1— 0), ^^C0(2 — 1), and ^^C0(1— 0) line emission in two selected arm 
regions in the central 6 kpc, a LVG analysis finds that the GMCs inside the spiral are similar 
to those observed in our Galaxy. The average kinetic temperature and H2 density (assuming 
virialized clouds) are T^in = 20K and n(H2) = 120 cm~^, respectively, at a linear resolution 
of 120 pc (290 pc for the southern region). Similarly, the derived conversion factor Xcn 



i s clos e to the Galactic value of 1.8 x 10^" cm ^K ^ km ^s recently derived by iDame et al. 



(I2OOII ). We interpret our results such that the physical properties of the CMC population do 



n ot change much as a function of environment in galaxy disks, consistent with the findings 



of iBolatto et al.l (120081 ) who analyzed CMC properties in local galaxies. 



Comparison of the derived molecular gas properties to the location of star formation 
shows no obvious trend suggesting that the massive star formation has no strong impact onto 
the GMCs inside the spiral arms. Comparison of the extinction inferred from the derived H2 
density and measured toward HII regions shows that the extinction toward the molecular 
clouds is at least a factor of 5 hig her with an average Ay of 15-30 mag than the average Ay 



of 3.1 inferred for HII regions by lScoville et al.l (I2OOII ) 
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Fig. 1. — Location of the two OVRO po intings (dark circle s) overlaid onto the integrated 
^^CO(l-O) line emission (contours) from lAalto et al.l ( 1l999l ) with short spacing correction 
(SSC) applied. The size of the circles corresponds to the primary beam at 3mm. The light 
gray rectangles outline the western (top) and southern region (bottom) that are studied. 
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Fig. 2. — The western region at 4.5" resolution in molecular line emission of (from left to 
right): i2C0(l - 0), i2C0(2 - 1), i3C0(l - 0), Ci^O (1-0), HCN (1-0), and HCO+(1-0). All 
maps shown are made from the SSC data cubes. The beam is shown in the bottom right 
corner of each panel. The angular offset is relative to the pointing center. 
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Fig. 3. — The western region at 2.9" resolution in molecular line emission of ^^C0(1 — 0) 
(left), ^^C0(2 — 1) {middle), and ^^C0(1 — 0) (right). All maps shown are made from the 
SSC data cubes. The beam is shown in the bottom right corner of each panel. The angular 
offset is relative to the pointing center. 
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Fig. 4. — The southern region at 4.5" resolution in molecular line emission of ^^C0(1 — 0) 
[top), ^^C0(2 — 1) {middle), and ^^C0(1 — 0) {bottom). All maps shown are made from the 
SSC data cubes. The beam is shown in the bottom right corner of each panel. The angular 
offset is relative to the pointing center. 
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Fig. 5. — Comparison of the location of the peak emission in the ^^C0(1 — 0) (gray contours) 
and ^^C0(1 — 0) (black contours) in the western arm at 2.9" resolution. The contours are 
in steps of 10% of the maximum value of 295Kkms~^ starting at 20% for the ^^C0(1 — 0) 
map and in steps of 10% of the maximum value of 37.3 K km s^^ starting at 30% for the 
^•^00(1 — 0) data. All maps shown are made from the SSC data cubes. 
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Fig. 6. — Comparison of the location of the peak emission in the ^^C0(1 — 0) (gray contours) 
and ^'^C0(1 — 0) (black contours) in the southern arm at 4.5" resolution. The contours are in 
steps of 10% of the maximum value of 201 K km s~^ starting at 30% for the ^^C0(1 — 0) map 
and in steps of 20% of the maximum value of 5.0Kkms~^ starting at 20% for the ^^C0(1 — 0) 
data. All maps shown are made from the SSC data cubes. 
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Fig. 7. — Moment maps of the ^^C0(2 — 1) line emission (at a resolution of 1.7" x 1.5" 
and short spacing corrected): integrated intensity map (left), velocity field {middle), and 
dispersion map {right). The contours in the intensity map are in steps of 10% of the peak 
value of 23.3 Jybeam~^kms~^, the iso-velocity contours start at 450km/s with steps of 
5 km/s, while the contours in the dispersion map start at 5 km/s with a step of 5 km/s. The 
beam is shown in the bottom right corner of each panel. The angular offset is relative to the 
pointing center. 
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Fig. 8. — Moment maps of the ^^C0(2 — 1) line emission (at a natural resolution of 1.7" xl.5" 
and short spacing corrected): intensity map (top), velocity field {middle), and dispersion 
map (bottom). The contours in the intensity map are in steps of 10% of the peak value of 
12.4 Jybeam~^kms~"'^ starting at 30%, the iso-velocity contours start at 490km/s with steps 
of lOkm/s, while the contours in the dispersion map start at 5km/s with a step of 5km/s. 
The beam is shown in the bottom right corner of each panel. The angular offset is relative 
to the pointing center. 
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Fig. 9. — Location of the positions studied in the western arm. Numbering of the positions 
is given in the left panel where the size of the circles corresponds to a resolution of 2.9". 
The derived line ratios R21 of the ^^C0(2 — 1) and ^^C0(1 — 0) emission and R13 of the 
^^C0(1 — 0) and ^^C0(1 — 0) emission for the individual positions is listed in the middle 
panel. Positions with empty circles have only upper limits in at least one line. The derived 
kinetic temperature Tkin and volume density of H2 n(if2) from the LVG analysis are given in 
the rightpanel. The integrated ^^C0(2 — 1) intensity (contours) is overlaid on the ^^C0(1— 0) 
intensity (gray-scale) . 
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Fig. 10. — Location of the positions studied in the western arm. Numbering of the positions 
is given in the top panel where the size of the circles corresponds to a resolution of 4.5". 
The derived line ratios R21 of the ^^C0(2 — 1) and ^^C0(1 — 0) emission and R13 of the 
^^C0(1 — 0) and ^^C0(1 — 0) emission for the individual positions is listed in the middle 
panel. Positions with empty circles have only upper limits in at least one line. The derived 
kinetic temperature Tkin and volume density of H2 n(iif2) from the LVG analysis are given 
in the bottom panel. The integrated ^^C0(2 — 1) intensity (contours) is overlaid on the 
^^C0(1 — 0) intensity (gray-scale). 
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Fig. 11. — Comparison of the ^^C0(2 — 1) line emission (same contours as Fig. [7]) to the 
stellar light distribution in the HST V band (left) as well as the Ha {middle left) and Paa 
{middle right) line emission mainly arising from HII regions and the MIPS 24yum continuum 
emission tracing warm dust heated by young stars {right). 
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Fig. 12. — Comparison of the ^^C0(2 — 1) line emission (same contours as Fig. [8]) to the 
stellar light distribution in the HST V band {top) as well as the Ha {middle top) and Paa 
{middle bottom) line emission mainly arising from HII regions and the MIPS 24/im continuum 
emission tracing warm dust heated by young stars {bottom). 
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Table 1. M 51a fields 



Pointing R.A. (J2000) Dec. (J2000) 



West 
South 



13:29:50.14 
13:29:55.17 



+47:11:25.27 
+47:10:47.41 



Note. — Pointing centers for two selected 
fields in the Ml spiral arm of M51a. 



Table 2. Observations and Set-up 



Pointing 


Transition 


Tracks 


Vlsr 
(kms-i) 


Av 
(kms-i) 


^ of Channels 


West 


i2CO(2 - 1) 


1.5L,1E,1H 


489 


2.6 


64 




i3CO(l-0) 


1C,1L,1.5H 


489 


2.7 


32 




Ci«O(l-0) 


1C,1L,1.5H 


488 


5.5 


16 




HCN(l-O) 


1C,2L,1.5H 


493 


6.8 


16 




HCO+(1-0) 


1C,2L,1.5H 


493 


6.7 


16 


South 


i2CO(2 - 1) 


1.5L,1E,1H 


517 


5.2 


32 




i3CO(l-0) 


1C,1L,1H 


517 


5.4 


16 




Ci8O(l-0) 


1C,1L,1H 


517 


5.5 


16 




HCN(l-O) 


0.5C,4L,2H 


523 


6.8 


16 




HCO+(1-0) 


0.5C,4L,2H 


523 


6.7 


16 



Note. — Here we list the final spectral resolution used and the vlsr of the 
central channel. The observed tracks are given as the number of full (~ 5hr 
on-source) tracks per configuration, e.g. 1.5L = 1.5 full tracks in L configuration. 
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Table 3. Parameters of the OVRO Data 



Pointing 


Line 


Weighting 


rms 
(mJybeam"^) 


beam 


West 


i2CO(2-l) 


na 


22 


1.67 X 1.63 




i2CO(2-l) 


ro 


26 


1.31 X 1.28 




i2CO(2-l) 


un 


27 


1.27 X 1.25 




^3CO(1-0) 


na 


17 


3.91 X 3.19 




i3C0(l - 0) 


ro 


19 


2.87 X 2.45 




Ciso(l-O) 


na 


13 


4.22 X 3.77 




HCN(l-O) 


na 


13 


4.76 X 3.98 




HCO+(1-0) 


na 


12 


4.67 X 3.93 


South 


i2CO(2-l) 


na 


19 


1.68 X 1.46 




i2CO(2-l) 


ro 


22 


1.33 X 1.21 




i2CO(2-l) 


un 


25 


1.28 X 1.20 




i3C0(l - 0) 


na 


17 


4.45 X 3.69 




Ci8O(l-0) 


na 


15 


4.63 X 4.18 




HCN(l-O) 


na 


12 


3.88 X 3.28 




HCO+(1-0) 


na 


11 


4.07 X 3.53 



Note. — Three different values are used for weighting of the uv data, 
namely a (Briggs) robust parameter of +5 (=na; 'natural weighting'), 
(=ro; 'robust weighting'), and -5 (=un; 'uniform weighting'). 
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Table 4. Sensitivity of the Short Spacing Corrected Data 

Pointing i2co(l-0) 12^0(2-1) i3co(l - 0) 01^0(1-0) 
(K) (K) (K) (K) 



West-2.9 


0.32 


0.13 


0.28 


— 


West-4.5 


0.23 


0.08 


0.09 


0.08 


South-4.5 


0.25 


0.06 


0.09 


0.09 



Note. — The rins noise in the short spacing corrected data cubes of 
the CO transitions. All CO data cubes have the same resolution of 2.9" 
(West-2.9) or 4.5" (West-4.5, South-4.5). 
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Table 5. Line Fluxes and Ratios at 2.9" resolution 



Position 


Offset 


i2C0(l - 


-0) 


1200(2 


-1) 


1300(1 


-0) 


^^50 


12C0(2-1) 
12CO(1-0) 


12CO(1-0) 
i3CO(l-0) 




(";") 


(K 


kms 


r') 


(K 


kms 


r') 


(K 


km; 


^~^) 


(kms-i) 






wl 


0;+15 






187 






116 






26 


60 


0.72 


5.5 


w2 


0;+ll 






150 






87 






21 


50 


0.60 


5.5 


w3 


-6;-6 






87 






103 






<19 


36 


1.10 


>9.5 


w4 


-2;+6 






193 






152 






39 


35 


0.65 


4.3 


w5 


+l;+5 






219 






135 






14 


60 


0.70 


7.5 


w6 


-l;+3 






204 






112 






22 


50 


0.60 


7.0 


w7 


+2;+2 






185 






137 






21 


46 


0.71 


4.5 


w8 


-1;0 






303 






157 






31 


40 


0.55 


7.5 


w9 


-2;-3 






222 






117 






<23 


45 


0.55 


>9.0 


wlO 


+5;-4 






102 






65 






<22 


43 


0.63 


>7.0 


wll 


+lr5 






185 






102 






31 


30 


0.45 


7.7 


wl2 


+3;-9 






128 






64 






23 


31 


0.45 


7.3 


wl3 


+4;-12 






142 






71 






28 


36 


0.47 


6.5 


wl4 


+6;-15 






215 






113 






47 


38 


0.42 


4.5 



Note. — Here we list the measured integrated line ffuxes for the different position in the western 
pointing. The lower case letters refer to the 2.9" resolution data that was used to derive the mea- 
surements. The X— and y— offsets are relative to the pointing centers given in Tab. [H The FWHM 
(vso) is derived from the 1^00(1 — 0) line. The line ratios were derived from scaling the 1^00(2 — 1), 
1^00(1 — 0), and 0^^0(1 — 0) spectra to the 1^00(1 — 0) spectrum. This method is more reliable 
than using the peak flux or the integrated flux in the case of low S/N data. Therefore the line ratios 
listed are not in perfect agreement with the ratio of the integrated line fluxes. 



Table 6. Line Fluxes and Ratios at 4.5" resolution 



Position 


Offset 


i2C0(l- 


-0) 


i2CO(2 


-1) 


i3C0(l- 


-0) 


Ci«0(l 


-0) 


^50 


12C0(2-1) 
12CO(1-0) 


12CO(1-0) 
i:^CO(l-0) 


12CO(1-0) 
Ci8O(l-0) 




(";") 


(K 


kms 


r') 


(K kms 


r') 


(K 


kms 


-') 


(K kms 


^-^) 


(kms-1) 








Wl 


0; + 15 






164 




91 






22 




<13 


64 


0.63 


5.1 


>15 


W2 


0;+ll 






151 




77 






19 




10 


55 


0.60 


5.8 


17 


W3 


-6;+7 






83 




89 






4 




<9 


45 


1.10 


10 


>io 


W4 


-2;+6 






184 




127 






25 




7 


39 


0.60 


6.0 


25 


W5 


+2;+4 






155 




115 






18 




<11 


55 


0.76 


5.7 


>14 


W6 


-1;0 






273 




136 






27 




6 


42 


0.53 


8.0 


25 


W7 


+l;-5 






170 




100 






18 




<7 


32 


0.50 


9.0 


>26 


W8 


+5;-10 






113 




67 






16 




<8 


39 


0.54 


7.8 


>14 


W9 


+7;-16 






178 




100 






25 




<8 


39 


0.45 


5.0 


>^ 


SI 


-14;+2 






44 




46 






<5 




<5 


26 


0.87 


>8.6 


^9 


S2 


-5;-3 






147 




57 






11 




<6 


34 


0.42 


11 


>23 


S3 


-2;+l 






134 




62 






7 




<7 


37 


0.48 


9.0 


>19 


S4 


+3;+l 






166 




59 






5 




<9 


46 


0.43 


11 


>19 


S5A 


+5;-l 






45 




26 






<6 




<6 


30 


0.52 


>7.5 


>8 


S5B 








86 




29 






<6 




<6 


30 


0.37 


>9.5 


>15 


S6A 


+13;0 






20 




24 






<6 




<8 


40 


0.90 


6.5 


>3 


S6B 








80 




36 






8 




<4 


22 


0.40 


7.3 


>20 


S7A 


+12;+5 






60 




26 






3 




<6 


33 


0.55 


10 


>10 


S7B 








32 




21 






2 




<5 


26 


0.71 


14 


>6 



Note. — Here we list the measured integrated line fluxes for the different position in the western and southern pointing. The 
capital letters refer to the 4.5" resolution data where all measurements were obtained. See caption of Tab. Ofor details. 
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Table 7. Results of the Large Velocity Gradient Analysis 



Position ^ 



dv 
dr 



J- kin 

(K) 



dr 



vir. 



N(H2) 
'1022cm-2l 



dr 



vir. 



nfHo 



cm 



^=c. 

dr 



Vir. 



(10 



20 



cm 



^=c. 

dr 



-2K-^km- 
vir. 



wl 


0.4 


22 


27 


1.5 


1.7 


340 


210 


0.8 


0.9 


w2 


0.4 


13 


13 


1.8 


2.1 


280 


160 


1.2 


1.4 


w4* 


0.5 


12 


13 


2.8 


3.5 


350 


230 


1.5 


1.8 


w5 


0.3 


35 


50 


1.2 


1.2 


250 


120 


0.6 


0.5 


w6 


0.3 


18 


22 


1.6 


2.4 


230 


120 


0.8 


1.2 


w7 


0.5 


16 


18 


2.1 


2.6 


380 


250 


1.1 


1.4 


w8 


0.3 


14 


16t 


2.9 


4.5 


200 


100 


0.9 


1.5 


wll 


0.3 


8 


lot 


3.2 


4.3 


180 


90 


1.7 


2.3 


wl2 


0.3 


8 


10 


2.2 


2.9 


190 


100 


1.7 


2.3 


wl4 


0.4 


5t 


5t 


10.9 


16.6 


280 


170 


5.1 


7.7 


S2 


0.2 


13 


15t 


1.4 


2.6 


120 


50 


0.9 


1.8 


S3* 


0.3 


13 


17 


1.3 


2.0 


160 


70 


1.0 


1.5 


S4* 


0.2 


13 


15 


1.5 


2.9 


120 


50 


0.9 


1.8 


S6B* 


0.3 


7t 


7t 


1.7 


3.0 


190 


100 


2.1 


3.7 


S7A 


0.2 


21 


24 


0.4 


0.7 


140 


60 


0.7 


1.2 


S7B^ 


0.3 


43 


50 


0.2 


0.3 


260 


110 


0.5 


0.8 


< value > 


0.3 


16 


20 


2.3 


3.3 


240 


120 


1.3 


2.0 



o 



'^A p2CO]/p3Qo] abundance of 60.0 was used. 
*Bright Ha emission present inside CO aperture. 
^Beani filling factor > 1. 

Note. — The kinetic temperatures Tkm and H2 column and volume densities N(H2) 
and n(H2) were derived from the LVG analysis described in sectionHl Typical uncertain- 
ties for the kinetic temperature are (75-100)% while the H2 densities have uncertainties 
of about 25% based on the (10-15)% calibration uncertainty for the CO hne data. A CO 
abundance of [12C0]/[H2] and [^^cOJ/i^^cQ] of 8.0 x 10"^ and 30.0, respectively, were 
assumed. The velocity gradient ^ was either kept fixed at 1.0 kms~^pc~^ or viriahzed 



gas was assumed where ^ = 3.1-v/ "^q4 . The corresponding conversion factors Xco 
from CO line temperature to H2 column density are also listed, the sta ndard Galactic 
conversion factor is 1.8 x 10^°cm~^K~^km~^s (e.g. iDame et al.ll200l[ l. The last row 
lists the average values of all regions analyzed. 
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